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ABSTRACT 

We report a direct and accurate measurement of the distance to the X-ray 
binary Cygnus X-1, which contains the first black hole to be discovered. The 
distance of 1.86l°;}? kpc was obtained from a trigonometric parallax measurement 
using the Very Long Baseline Array. The position measurements are also sensitive 
to the 5.6 d binary orbit and we determine the orbit to be clockwise on the 
sky. We also measured the proper motion of Cygnus X-1 which, when coupled 
to the distance and Doppler shift, gives the three-dimensional space motion of 
the system. When corrected for differential Galactic rotation, the non-circular 
(peculiar) motion of the binary is only about 21 km s~^, indicating that the 
binary did not experience a large "kick" at formation. 

Subject headings: astrometry — black hole physics — stars: distances, individual 
(Cygnus X-1) — X-rays: binaries 



1. Introduction 

In this paper (Paper I), we resolve the long-standing problem of the distance to Cygnus X- 
1 by measuring its trigonometric parallax, ie, triangulating using the Earth's orbit as one 
leg of the triangle and measuring the change in its apparent position as the Earth orbits 
the Sun. Our results are based on observations of compact radio emission with the National 
Radio Astronomy Observatory'^ Very Long Baseline Array (VLBA). Previously, only one 
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trigonometric parallax has been published for a binary containin g a black hole: V404 Cyg 



an X-ray binary with a low-mass (M < IMq) companion star (IMiller- Jones et al. 1 120091 ) . 



Our result is the first trigonometric parallax for an X-ray binary with a massive companion 
(M > IOMq); an accurate distance is key to achieving secure and detailed models of the 
physical properties of the binary system. 

Cygnus X-1 was the first black ho l e (BH) candid ate to be established via dynamical 
observations ( IWebster fc Murdin I Il972t iBolton 1 119721 ) . It is one of the most well-studied 
black hole systems with excellent optical and X-ray data. However, modeling this system is 
limited by its very uncertain dist ance. Distance es timates span a very large range, with most 
values between 1.8 and 2.4 kpc ( iZiolkowski II2005I ). The Cygnus X-1 distance estimate with 



the sm allest formal uncertainty is 2.14 ±0.07 kpc by lMassey. Johnson &: Degioia-Eastwood 



( 119951 ) and is based on photometry and spectroscopy of NGC 6871 (an open cluster in the 
Cyg 0B3 complex). While this distance estimate is very precise, its true accuracy (including 
systematic sources of error involving extinction, metallicity corrections, and stellar modeling) 
is perhaps ±15%. Also, it (and many other distance estimates) relies on the association of 
Cygnus X-1 with the Cyg 0B3 complex. While such an association is supp orted by similar 



stellar (Hipparcos) and Cygnus X -1 (VLBI) proper motion measurements (ILestrade et al. 



19991 : iMirabel &: Rodrigues II2003I ). it still must be considered tentative as this line-of-sight 
through the Milky Way is very "crowded," passing nearly parallel to the Local spiral arm, 
which extends at least sev eral kpc from the Sun, and then through the Perseus and Outer 
arms at greater distances (IReid et al. Il2009al ). 



2. Observations and Data Analysis 

Our observations were conducted under VLBA program BR141. We observed Cygnus X- 
1 and two background sources over 10-hour tracks at five epochs: 2009 January 23, April 
13, July 13, October 31 and 2010 January 25. These dates well sample the peaks of the 
sinusoidal trigonometric parallax signature in both Right Ascension and Declination. This 
sampling provides near maximum sensitivity for parallax detection and ensures that we can 
separate the secular proper motion (caused by projections of Galactic rotation, as well as 
any peculiar motion of Cygnus X-1 and the Sun) from the sinusoidal parallax effect. Tabled] 
lists the positions of Cygnus X-1 and the background continuum sources, as well as their 
angular separations from Cygnus X-1. 



Generally, data calibration followed similar procedure s as for parallax obs ervations of 
continuum sources in the Orion nebular cluster at 8.4 GHz (IMenten et al. 1120071 ). We placed 
observations of well-known strong sources near the beginning, middle and end of the obser- 
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vations in order to monitor delay and electronic phase differences among the intermediate 
frequency bands. In practice, however, we found minimal drifts and used only a single scan 
of J2005+7752 for this calibration. 

In order to calibrate atmospheric propag ation delays, we scheduled "geodetic blocks" 



consisting of typically a dozen ICRF sources (IFey et al.l 120041 ) obser ved over a wide range 
of source azimuths and elevations as describe in iReid et al. I ( l2009bl ). We observed in left 
circular polarization and spread eight 16-MHz bands so as to span ^ 500 MHz of bandwidth. 
We corrected the interferometer data for ionospheric delays using glo bal models of the total 



electron content obtained from global positioning system observations ( jWalker fc Chatter jee 



20001 ). Residual multi-band delays and fringe rates were measured and used as data to solve 
for clock offsets and drift rates and unmodeled zenith atmospheric propagation path delays. 
Later these were removed from the rapid-switching (parallax) data. 

The rapid-switching observations employed four adjacent frequency bands of 16-MHz 
bandwidth and recorded both right and left circularly polarized signals. The four (dual- 
polarized) bands were centered at frequencies of 8.425, 8.441, 8.457 and 8.473 GHz. We 
alternated between two ~ 11 min blocks, each consisting of observations of Cygnus X- 
1 and one of the background sources. Within a block, we switched sources every 40 s, 
typically achieving 30 s of on-source data. We used the background sources J1953-I-3537 and 
J1957+3338 for phase-referencing. 

We calibrated the correlated data using the NRAO Astronomical Image Processing Sys- 
tem (AIPS). The calibration sequence included four steps. The first step involved correction 
of interferometer delays and phases for the effects of diurnal feed rotation (parallactic angle), 
for errors in the values of the Earth's orientation parameters used at the time of correlation, 
and for any small position shifts in the a priori coordinates of either the maser or back- 
ground sources. Since the VLBA correlator model includes no ionospheric delays, we used 
global total electron content models to remove ionospheric effects. At this point, we also 
corrected the data for residual zenith atmospheric delays and clock drifts (determined from 
the geodetic block data). 

In the second step, we adjusted interferometer visibihty amplitudes for the small (few %) 
effects of biases in the threshold levels of the data samplers at each station. We also entered 
system temperature and antenna gain curve information into calibration tables, which allows 
conversion of correlation coefficients to flux density units. In the third step, we performed a 
"manual phase-calibration" to remove delay and phase differences among all eight 16-MHz 
bands. This was done by selecting one scan on a strong calibrator, fitting fringe patterns 
to the data for each frequency band, and then shifting delays and phases to remove offsets. 
The fourth calibration step involved determining antenna-based complex gains from each 
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of the two background sources and interpolating and removing those from the Cygnus X-1 
observations. For most antennas at most times the phases were easily "connected." However, 
when interferometer phase differences between adjacent background source scans exceeded 
60°, the data between those times were discarded. 

We then imaged the calibrated data with the AIPS task IMAGR. The synthesized 
"dirty" beam typically had a full- width at half-maximum of 1.55 x 0.86 mas elongated 
approximately in the north-south direction. We adopted a round convolving-beam of 1.0 
mas to construct CLEAN images. Example images from the middle epoch are shown in 
Fig. [U Positions were determined by fitting elliptical Gaussian brightness distributions to 
the images. For Cygnus X-1, which shows a core-jet structure, we restricted the fitted pixels 
to those close to the bright core, typically within 1 mas of the peak emission. The peak 
brightness of Cygnus X-1 ranged between 4 and 9 mJy beam~^ among our observations. 



3. Parallax and Proper Motion Fitting 

The change in position of Cygnus X-1 relative to a background continuum source was 
then modeled by the parallax sinusoid in both coordinates, completely determined by one 
parameter (the parallax), and a secular proper motion in each coordinate (see Fig. The 
model included the effects of the ellipticity of the Earth's orbit. The weighting of the 
data in the parallax and proper motion fitting is complicated because the formal position 
uncertainties are often unrealistically small, since a priori unknown sources of systematic 
error often dominate over random noise. The north-south components of relative positions 
often have greater uncertainty than the east-west components, because the interferometer 
beams are generally larger in the north-south direction and systematic errors from unmod- 
eled atmospheric delays usually are more strongly correlated with north-south positions 



flHonma. Tamura fc Reid \\20Q& \. 



In order to allow for, and estimate the magnitude of, systematic errors, we assigned 
independent "error floors" to the east-west and north-south position data and added these 
floors in quadrature with the formal position-fitting uncertainties. Trial parallax and proper 
motion fits were conducted and a reduced xl (per degree of freedom) statistic was calculated 
separately for the east-west and north-south residuals. The error floors were then adjusted 
iteratively so as to make f» 1.0 for each coordinate. This procedure resulted in error floors 
of 0.08 and 0.16 mas for the eastward and northward position measurements, respectively. 
The magnitudes of these error floors are reasonably consistent with those obtained for other 



parallax targets observed at 8.4 GHz, e.g. iMenten et al. I (120071 ). and are probably dominated 



by unmodeled ionospheric delays. Any component of variability in the centroid of the core-jet 
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Table 1. Source Positions 



Source 


R.A. (J2000) 


Dec. (J2000) 


^sep 


P.A. 




(h m s) 


(d ' ") 


(deg) 


(deg) 


Cygnus X-1 


19 58 21.672762 


35 12 05.72511 






J1953+3537 ... 


19 53 30.875712 


35 37 59.35927 


1.1 


-66 


J1957+3338 ... 


19 57 40.549923 


33 38 27.94339 


1.6 


-175 



Note. — The position for Cygnus X-1 at 8.45 GHz comes from the par- 
allax fitting at epoch 2009.570. This position is relative to the unweighted 
average of the two background sources, which gave coordinate differences 
of about ± 0.25 mas, consis tent with their ICRF catalog uncertainties 
of 0.5 mas (IFey et al.ll2004j ). Angular offsets (6'scp) and position angles 
(P.A.) east of north relative to Cygnus X-1 are indicated in columns 4 
and 5. 




Fig. 1. — Images of Cygnus X-1 and the two background continuum sources from the middle 
epoch (2009 July 13). Source names are in the upper left corner and restoring beams are in 
the lower left corner of each panel. All contour levels are integer multiples of 1 mJy beam~^ 
for Cygnus X-1 and 15 mJy beam~^ for the background sources. Zero contours are suppressed 
for clarity; negative brightness levels were below the first contour level for all images. 
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position of Cygnus X-1 caused by changing jet opacity must be less than ^ 0.1 mas. 

The fitting results are presented in Table [21 First, we attempted separate parallax 
fits for each background source measurement. The parallaxes for Cygnus X-1 for the two 
background sources differ by only 0.071 mas and are consistent with their individual formal 
uncertainties of about 0.06 mas. This agreement suggests that neither atmospheric mis- 
modeling nor potential unresolved structural variations in the background sources contribute 
position errors greater than the error floors. 

Next, because of the agreement of the Cygnus X-1 parallaxes and proper motions for 
the two background sources, we did a combined fit with common parallax and proper motion 
parameters for both sets of data. (In the combined fit, we did allow for different position 
offsets for Cygnus X-1 with respect to each background source, since the a priori positions 
for these sources were not known to better than ^ 0.5 mas.) 

Finally, we investigated the effect of the Cygnus X-1 binary orbit on the parallax fits. All 
orbital pa rameters are known with reasonable accuracy and were consid ered as independent 



variables (IBrocksopp et al. Ill999l : iGies et al. Il2008t lOrosz et al. Il201ll ). with the exception 



of the direction of the orbital rotation on the sky. The binary orbital model we adopted 
had an orbital period of 5.599829 d, superior conjunction at JD 2441874.707 d, a radius of 
the radio source about the center of mass of 0.128 AU, an inclination of 36° (ie, the angle 
between the orbit angular momentum vector and our line-of-sight). (Using the revised value 
of i = 27°, determined in Paper II, results in no significant difference in parallax and proper 
motion fits.) We assumed an orbital angular momentum vector paral lel to the radio jet in 



Cygnus X-1, observed in our images (see Fig. [T]) and previous studies (jStirling et al. 1120011 ) 
to be about —26° east of north. This leaves two possible orbital directions, clockwise (CW) 
and counterclockwise (CCW) on the sky, to be determined by the data. 

Parallax and proper motion fits for both orbital directions are listed in Table [2l The 
value for the CW orbit was 8.49 (for 13 degrees of freedom), which was significantly lower 
than the value of 16.81 for the CCW orbit. The agreement of the data with the CW orbit 
model can be seen in the position data in Fig. [31 especially in the more accurate eastward 
offset data. Thus, our astrometric observations appear to be sensitive to the position changes 
caused by orbit of the radio source (presumably the black hole) about the center of mass of 
the binary system, and we have directly determined the direction of the orbit to be clock-wise 
on the sky. 

Adopting the CW orbit, we re- fitted the data allowing the orbital radius to vary and 
found a radius of 0.18 ± 0.09 AU. Thus, the current VLBA data independently gives a 2a 
estimate of the radius of the compact radio source about the center of mass of the binary. 



- 7- 




2 0-2 2009 2010 2009 2010 

East Offset (mas) Epoch (years) Epoch (years) 



Fig. 2. — Parallax and proper motion data and fits. Plotted are position offsets of Cygnus X- 
1 relative to the two background sources: J1953+3537 (red triangles) and J1957+3338 (blue 
squares ). The degree-scale separations of the absolute positions given in Table [1] have been 
removed. The yearly sinusoidal position variations (the parallax signature) caused apparent 
position shifts of Cygnus X-1 as the Earth orbits the Sun. Left Panel: Positions on the sky 
with first and last epochs labeled. Data for the two background source are offset horizontally 
for clarity. The expected positions from the parallax and proper motion fit are indicated 
(circles). Middle Panel: east (top lines) and North (bottom lines) position offsets and best fit 
parallax and proper motions fits versus time. Data for the two background sources and two 
coordinates are offset vertically for clarity. Right Panel: Same as the middle panel, except 
the with the best fit proper motions have been removed, allowing all data to be overlaid 
and the effects of only the parallax seen. The north offset data have been offset vertically 
(below) the east offset data for clarity. 
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Table 2. Cygnus X-1 Parallax & Proper Motion Fits 



Cygnus X-1 Orbit 


Background 


Parallax 








N 




Source 


(mas) 


(mas y~^) 


(mas y~"^) 






None 


J1953+3537 


0.582 ± 0.055 


-3.74 ±0.10 


-6.56 ±0.19 


4.57 


5 




J1957+3338 


0.511 ±0.061 


-3.66 ±0.11 


-6.28 ±0.21 


5.54 


5 




Combined 


0.547 ±0.041 


-3.70 ± 0.08 


-6.42 ±0.14 


12.82 


13 


CCW 


Combined 


0.528 ±0.047 


-3.68 ±0.09 


-6.36 ±0.17 


16.81 


12 


cw 


Combined 


0.539 ±0.033 


-3.78 ± 0.06 


-6.40 ±0.12 


8.49 


12 



Note. — Combined fits used a common parallax and proper motion parameters for both 
background sources. The proper motion components are defined as /ix — /J'a cos 5 and /ly — fj,s, 
where a is Right Ascension and 5 is Declination. The binary orbital model assumed an orbital 
period of 5.599829 d, superior conjunction at JD 2441874.71 d, a radius of the radio source 

about the center of mass of 0.071 mas, an inclination of 36°, and spin-axis of —26° east of 
north. We tried two orbital direction: clockwise (CW) and counter-clockwise (CCW) on 
the sky. Columns 6 and 7 list the post-fit values and the number of degrees of freedom 
(N). Parameter uncertainties have been re-scaled to those appropriate for a xl (p^r degree 
of freedom) of unity. 
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With more observations in the future, we could refine the values of some of the Cygnus X-1 
orbital parameters. 



We adopt a Cygnus X-1 parallax of 0.539 ± 0.033 mas from the fit with the lowest xt 
value. That fit combined the data from both background sources and used the orbital model 
with clockwise projected motion on the sky. This parallax for Cygnus X-1 corresponds to 
a distance of 1.86^q;}^ kpc. At this distance, the proper motion of Cygnus X-1 corresponds 
to —33 km s~^ and —56 km s~^ eastward and northward, respectively. Completing the full 
space velocity, th e average heliocentric radial velocity of Cygnus X-1 is about —5 km s~^ 
flCies et al. Iboosi ). 



4. Discussion 



Since we have measured the position, parallax, proper motion, and others have measured 
the line-of-sight velocity, we have full three-dimensional position and velocity information for 
Cygnus X-1. We can convert this velocity vector from the equatorial heliocentric reference 
frame in which they are measured into a Galactic reference frame. A convenient Galactic 
frame is one rotating with a circular velocity B(-R) at the position of the source: ie, a "local 
st andard of rest" at t he location of Cygnus X-1. This transformation is described in detail 
2009ah . 



m 



Reid et al. 



Conversion to a rotating Galactic reference frame depends somewhat upon the solar 
motion and Galactic parameters. We adopt the distance of the Sun from t he G alactic 
center of Rq = 8.25 kpc as a rounded average of the valu es q uoted bv iReid I (119931) for a 
variety of "classical" techniques and by iGhez et al. I fcOOSl ) and iGillessen et al I (120091 ) from 
the orbits of stars about the supermassive black hole at the Galactic center. The rotation 
speed of the local standard of rest is Qq ~ 240 km s~^, which is based on this distance and 
the proper motion of Sgr A* (the super- massive blac k hole at the Galactic cen ter), which 
directly measures (6o + Vq)/Rq = 30.2 km s~^ kpc~^ jReid fc Brunthaler" 2004). Us i ng th e 
newly proposed Solar motion of V© ~ 11 km s~^ by ISchonrich. Binney fc Dehnen I ( 120101 ). 



and assuming a locally fiat rotation curve, i.e. Q{R) = Go, we find a "peculiar motion" 
(relative to circular motion) for Cygnus X-1 of 17 ± 1 km s~^ toward the Galactic center, 
— 11 ± 3 km s~^ in the direction of Galactic rotation, and 6 ± 1 km s^^ toward the North 
Galactic Pole. Thus, the total magnitude of the peculiar motion of Cygnus X-1 is only 
21 km s~^, indicating that the binary system did not experience a large "kick" when the 
black hole formed. This result s upports the conjecture that the black ho le may have formed 
without a supernova explosion (IKatz Ill975l : iMirabel fc Rodrigues II2003I ). 
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2009 2009.25 2009.5 2009.75 2010 



Epoch (years) 

Fig. 3. — Parallax fit including the effects of a clock-wise (on the sky) binary orbit for 
the radio source relative to the center of mass of Cygnus X-1. East (top) and north (bot- 
tom) position offsets and best fitting model (lines) versus time. Plotted are the positions 
of Cygnus X-1 relative to the two background sources: J1953+3537 (red triangles) and 
J1957+3338 (blue squares). Clearly seen is the annual parallax sinusoidal signature (with 
the proper motion removed) and the "high frequency" oscillations caused by the 5.6 day or- 
bit of the radio source (the black hole) about center of mass of the binary system. Blow-ups 
of the eastward (top small boxes) and northward (bottom small boxes) offsets each cover 0.1 
year in time and 0.25 mas in offset. 
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Recently, IXiang et al. I ( 120111 ) modeled the scattering by interstellar dust grains of x- 
rays from Cygnus X-1. They considered 18 models describing the characteristics of dust and 
selected six models with iron column densities consistent with a value of 1.62 x 10^^ cm~^ 
(Lee et al., in preparation). Modeling the x-ray light-curves, three of the six yielded distances 
between 3.15 and 3.59 kpc, and these models were discarded as being incompatible with our 
trigonometric parallax. The remaining three models yielded distances between 1.72 and 1.90 
kpc, compatible with our parallax. Thus, having an accurate trigonometic parallax can be 
important for the study of x-ray scattering from time varying sources. 

The black hole in Cygnus X-1 is fed gas via the stellar wind emanating from its 0-type 
companion, which underfills its Roche equipotential lobe. The degree to which the star 
underfills its Roche lobe must be tightly constrained in order to obtain a useful dynamical 
mod el of the system, an d accurate knowledge of the source distance is necessary to accomplish 
this (lOrosz et al. 1120071 ). Using the accurate distance, reported here, we have constructed 
such a model and obtained precise values for the mass, M, of the black hol e, the orbital 
inclination angle, i, and other system parameters, as reported in Paper II (jOrosz et al. 



20111). With the req uired, accurate values of M, i and D in hand, we proceed in Paper 
III (jCou et al 2011 ) to model the X-ray continuum spectrum and determine the spin of 
the black hole. BH spin can be determined because the measured X-ray spectrum from the 
accretion disk is dominated by radiation coming from very near the BH. The temperature and 
luminosity of the X-ray emission increase greatly with increasing BH spin, as the inner radius 
of the accretion disk, R^^, shrinks toward the BH event horizon. R^^ can be determined by 
fitting the X-ray continuum emission from the BH's accretion disk to a fully relativis tic model 



if the black hole mass, distance and accretion disk inclination angle are known (ILi et al. 



20051 ). Because an astronomical black hole (which is believed to have negligible net electric 
charge) is completely described by specifying its mass and spin, our measurements of these 
quantities provide a complete description of the black hole in Cygnus X-1. 



Facilities: VLBA 
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